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ABSTRACT 

A  theory  of  interval  iteration,  based  on  a  few  simple  assumptions,  is 
given  for  the  fixed  point  problem  for  operators  in  partially  ordered 
topological  spaces.  A  comparison  of  interval  with  ordinary  iteration  is  made 
which  shows  that  their  properties  are  converse  in  a  certain  sense  with  respect 
to  existence  or  nonexistence  of  fixed  points.  The  theory  of  interval 
iteration  is  shown  to  hold  without  modification  if  the  computation  is 
restricted  to  a  finite  set  of  points,  as  in  actual  practice.  In  this  latter 
case,  interval  iteration  is  shown  to  converge  or  diverge  in  a  finite  number  of 
steps,  for  which  an  upper  bound  is  given.  By  the  introduction  of  a  suitable 
iteration  operator,  the  method  of  interval  iteration  is  extended  to  the 
problem  of  solution  of  equations  in  linear  space-s. 
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SIGNIFICANCE  AND  EXPLANATION 

Interval  iteration  is  a  powerful  computational  method  with  applications 
to  the  solution  of  linear  and  nonlinear  systems  of  equations,  integral 
equations  (see  MRC  TSR  #2128),  and  differential  equaitons.  Because  of  this 
diversity,  it  is  helpful  to  have  a  general  theory  to  act  as  a  guide  to 
specific  applications.  This  report  gives  such  a  general  theory,  based  on  only 
a  few  simple  assumptions  to  give  it  as  wide  a  scope  as  possible.  The  notions 
of  convergence  and  divergence  of  interval  iteration  are  defined  and  their 
implications  with  regard  to  existence  or  nonexistence  of  solutions  are 
presented.  It  turns  out  that  convergence  of  interval  iteration  is  assured  if 
the  initial  interval  contains  a  solution,  and  the  iteration  in  this  case  can 
provide  improved  lower  and  upper  bounds  for  the  solution.  Divergence  of 
interval  iteration,  on  the  other  hand,  shows  that  the  initial  interval  did  not 
contain  a  solution,  and  thus  may  be  discarded.  These  properties  are  compared 
with  the  well-known  results  for  ordinary  iteration,  and  shown  to  be  converse 
in  a  certain  sense.  Remarks  on  how  to  use  interval  iteration  in  applications 
are  made,  and  the  toeory  is  also  shown  to  hold  if  the  computation  is 
restricted  to  a  finite  set  of  elements,  as  on  an  ordinary  computer. 
Construction  of  interval  extensions  which  are  exactly  representable  in  the 
finite  case  is  explained,  and  it  is  observed  that  interval  iteration  will 
converge  or  diverge  in  a  finite  number  of  steps  in  actual  computation,  and  an 
upper  bound  is  obtained  for  this  number  which  depends  only  on  the  initial 
interval  chosen.  Now  that  microprogrammed  interval  operations  can  be 
performed  with  speeds  comoarable  to  floating  point  arithmetic,  interval 
iteration  can  be  applied  to  a  wider  class  of  problems  to  obtain  solutions  of 
known  reliability. 

The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


A  THEORY  OF  INTERVAL  ITERATION 


L.  B.  Rail 

1.  A  fixed  point  problem.  Suppose  that  $  is  an  operator  which  maps  a  partially 
ordered  topological  space  S  into  itself.  A  fixed  point  y*  £  S  of  $  satisfies  the 
equation 

(1.1)  y  =  4  (y) 

solving  this  equation  is  called  the  fixed  point  problem  for  4>  in  S. 

The  partial  ordering  relation  in  S  will  be  denoted,  as  usual,  by  Elements 

y_,  y  of  S  such  that  ji  ;  y  define  an  interval  Y  =  [^  ,  y]  in  S,  which  is  the  nonempty 
set 

(1.2)  Y  =  (£  ,  y]  =  (y  |  ^  s  y  s  y,  y  e  s}. 

The  elements  y  are  called  respectively  the  lower  and  upper  endpoints  of  Y.  The 
set  of  all  intervals  in  S  is  denoted  by  IS.  The  elements  y  of  S  are  identified  with 
the  corresponding  degenerate  intervals  y  =  (y  ,  y]  which  have  equal  endpoints. 

The  following  assumptions  will  be  made  concerning  the  topology  of  S: 

(i)  Intervals  are  closed  subsets  of  S; 

(ii)  Each  nondegenerate  interval  contains  a  limit  point  of  countable  order. 
The  last  assumption  means  that  if  Y  is  nondegenerate,  then  it  contains  a  limit 
point  i  such  that  each  neighborhood  of  l  contains  at  least  a  countable  number  of 
points  of  Y  different  from  i.  This  technical  property  is  called  N^-compactness  by 
Sierpirfski  [10) . 


and  (ii)  4  is  Inclusion  monotone,  that  is, 

(2.2)  YCz  •  WICJBI,  Y.Z  £  IS. 

♦  is  called  an  interval  extension  of  an  interval  operator  f  if  F(Y)  C  4 (Y)  for 
Y  e  is  and  (2.2)  holds. 

Definition  2.1.  Given  an  initial  interval  Yp  and  an  interval  operator  4,  the 

sequence  (Y  }  defined  by 
n 

(2.3)  Y  ^  -  Y  n*(Y),  n  -  0,1,2,..., 

n+i.  n  n 

is  said  to  be  generated  by  interval  Iteration  starting  from  Y^ . 

Note  that,  from  (2.3), 

(2.4)  y(|:y1dyj3...  i 

hence,  interval  iteration  generates  a  nested  (or  descending)  sequence  of  closed  sets. 
Definition  2.2.  The  interval  iteration  process  (2.3)  is  said  to  diverge  if 

(2.5)  Yn  ■  0  (empty) 

for  same  positive  integer  Nj  otherwise, 

« 

(2.6)  Y*  -  n  Y 

A  n 
n»0 

is  nonempty  by  the  Cantor  theorem  ([10],  pp.  34-35),  and  the  interval  iteration  is 
said  to  converge  to  the  limit  Y*  «  {Y^ }  given  by  (2.6) . 

Thus,  according  to  this  definition,  convergence  and  divergence  of  interval 
iteration  have  a  converse  relationship,  as  one  would  expect. 

Theorem  2.1.  If  4  is  an  interval  extension  of  the  operator  4  in  S,  and  the 
initial  interval  YQ  contains  a  fixed  point  y*  of  4,  then  the  interval  iteration  (2.3) 
converges  j  furthermore , 

m 

(2.7)  y»  €  Y»  -  lim  {Y  )  »  n  y  . 

n-*®  n  *.  n 

n-0 

Proof.  This  follows,  as  is  well  known  [3],  [6],  from  (2.3)  and  mathematical 
induction,  since  y*  6  Y^  “  4(y*)  -  y*  S  ifY^) ,  and  thus  y*  €  Y^+1  -  Y^fttJY^)  ?  0.  (JED 
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The  contrapositive  of  Theorem  2.1  is  the  following  result,  which  is  more  in¬ 
cisive  . 

Theorem  2.2.  If  4  is  an  interval  extension  of  4,  and  the  interval  iteration 
(2.3)  diverges,  then  the  initial  interval  Y^  contains  no  fixed  points  y*  c'  4. 

This  assertion  was  noted  by  Nickel  f8)  in  connection  with  an  interval  version 
of  Newton's  method. 

Observe  that  convergence  of  interval  iteration  is  a  necessary,  not  sufficient, 
condition  for  the  existence  of  a  fixed  point  y*  £  Y^  of  4;  divergence,  on  the  other 
hand,  is  a  sufficient  condition  for  nonexistence  of  fixed  points  of  4  in  YQ. 

3.  Comparison  with  ordinary  iteration.  The  ordinary  iteration  method 

(3.1)  yn+1  =  4<yn)  .  n  ■=  0,1,2,..., 

is  often  used  to  attemot  to  generate  a  sequence  {y  )  which  converges  to  a  fixed  point 
*■  n 

y*  of  4,  starting  from  some  initial  point  y  .  If  4  is  continuous  in  the  topology  for 

S,  in  which  also  lijn{y  }  »  y*  £  S,  then  y*  will  be  a  fixed  point  of  4,  and  if  Y„  is 
n-*“  11  0 

a  closed  subset  of  S  such  that  {y  }  C  Y,, ,  then  y*  £  Y„.  On  the  other  hand,  if  Y  is 

n  0  0  0 

a  subset  of  S  which  does  not  contain  a  fixed  point  of  the  continuous  operator  4» 
then  the  sequence  generated  by  the  iteration  (3.1)  cannot  converge  to  a  point  of  Yq. 
The  first  alternative  will  be  called  convergence  into  Y^,  and  the  second  divergence 
from  Yg.  On  the  basis  of  these  definitions  tod  the  corresponding  concepts  for  inter¬ 
val  iteration  given  in  Definition  2.2,  a  comparison  of  ordinary  and  interval  iteration 
is  shown  in  Figure  3-1. 


Ordinary  Iteration  y  _  =  4(y  ) . 

n+ 1  n 

Interval  Iteration  Y  _  *  Y  n  <t>Y  . 

n+I  n  n 

Convergence  (into  YQ)  * 

Existence  (y*  E  YQ) 

Existence  (y*  £  YQ)  “ 

OO 

Convergence  (Y  =  n  y  yi  0) 

_  n 
n»0 

Nonexistence  (y*  £  Y  )  * 

Divergence  (from  Y  ) 

Divergence  rsome  Y  *  0)  ■* 

N 

Nonexistence  (y*  £  Y^) 

Figure  3.1.  Ordinary  and  interval  Iteration  Compared. 


Thus,  interval  iteration  stands  in  a  converse  relationship  to  ordinary  itera¬ 
tion  under  the  above  assumptions.  It  is  worth  noting  that  in  metric  spaces  S,  the 
convergence  of  the  ordinary  iteration  process  (3.1)  often  depends  on  being  able  to 
choose  the  initial  point  yQ  "close"  to  the  fixed  point  y*,  and  that  some  operators  ♦ 
have  fixed  points  y*  which  repel  the  iteration  sequence  (yn>  for  all  yQ  f  y* .  The 
convergence  of  interval  iteration,  on  the  other  hand,  follows  if  the  initial  interval 
Yq  is  "large  enough"  to  contain  a  fixed  point  y*  of  $ . 

4.  Applications  of  interval  iteration.  Interval  iteration  may  be  applied  in  sev¬ 
eral  ways  to  the  fixed  point  problem  (1.1). 

1*.  Suppose  that  the  interval  YQ  is  known  to  contain  a  fixed  point  y*  of  <f, 
perhaps  on  the  basis  cf  a  nonconstructive  fixed  point  theorem.  In  this  case,  the 
interval  iteration  (2.3)  will  converge,  and  may  be  used  to  obtain  lower  and  upper 
bounds  for  y*,  namely, 

(4.1)  s  y*  s  yn,  n  «  0,1,2,..., 

where  Y  «  [»  ,  y  ] ,  and  in  the  limit , 

n  n 

(4.2)  s  y*  s  y*, 

where  Y»  -  [y*  ,  y»]  is  the  limit  (2.6)  of  {Y  }. 

n 

The  bounds  (4.1)  give  improved  results  as  long  as  the  inclusions  (2.4)  are 

strict.  However,  if  Y„.,  -  Y„  for  some  positive  integer  N,  then 
N+l  N 

m 

(4.3)  Y»  -  n  Y  -  Y  ; 

A  n  n 
n«0 

this  is  called  finite  convergence  of  the  interval  iteration  (2.3).  If  finite  con¬ 
vergence  takes  place,  then  the  best  lower  and  upper  bounds  obtainable  for  y*,  start¬ 
ing  from  Yq  -  .  yQ) .  are 

(4.4)  ij  «  S  y*  S  yH  -  y*  . 

2*.  If  it  is  not  known  whether  or  not  YQ  contains  a  fixed  point  y*  of  if,  in¬ 
terval  iteration  r.ay  still  be  useful  in  one  of  the  following  ways. 
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(i)  As  long  as  the  interval  iteration  is  producing  intervals  which  de¬ 
crease  at  each  step  (strict  inclusion  holds  in  (2.4)),  then  an  existence  test  that 
fails  because  Y  ,  Y^  Y^  are  "too  large"  may  succeed  for  Y^.  This  Y^  may  then 

be  taken  as  the  initial  interval  YQ,  and  one  has  the  favorable  case  1"  discussed  above 

(ii)  If  the  iteration  produces  an  empty  intersection  (divergence) ,  then 
this  establishes  conclusively  that  the  initial  interval  Y^  does  not  contain  a  fixed 
point  y*  of  $,  so  that  this  interval  may  be  excluded  from  further  consideration. 

There  is,  of  course,  a  third  possibility: 

(iii)  The  interval  iteration  (2.3)  leads  only  to  an  interval  Y°  in 
which  no  conclusive  assertion  about  existence  or  nonexistence  of  a  fixed  point  y* 
of  <p  is  available.  (Y°  may  be  the  limit  Y*  if  obtained  in  a  finite  number  of  steps, 
or  otherwise.) 

Possible  alternatives  in  this  situation  include  partition  of  the  resulting  in¬ 
terval  Y°  into  subintervals  for  further  examination,  a  strategy  developed  by  Moore 
and  Jones  [7]  ,  or  accsptance  of  Y°  as  a  "generalized"  or  "pseudosolution’*  (relative 
to  the  initial  interval  Yq)  of  the  fixed  point  problem.  This  latter  choice  may  be 
useful  in  the  development  of  an  interval  version  of  regularization  of  solutions  of 
ill-posed  problems. 

5.  Interval  iteration  on  a  grid.  The  ordinary  iteration  process  (3.1)  is  a  poor 
model  of  what  actually  occurs  in  computation,  since  it  is  usually  impossible  to  car¬ 
ry  out  the  indicated  transformations  exactly.  Interval  iteration,  on  the  other  hand, 
is  readily  adaptable  to  actual  machine  computation,  and  its  theory  can  be  preserved 
intact . 

Suppose  that  G  (called  a  grid)  is  a  finite  subset  of  the  space  S.  Here,  one 
may  think  of  the  set  of  numbers  which  have  exact  representations  on  a  given  computer, 
and  finite  Cartesian  products  of  such  a  set.  The  subset  of  IS  consisting  of  intervals 
with  endpoints  in  G  will  be  denoted  by  IG,  that  is, 

(5.1)  IG  =  {[a  ,  b]  j  a,b  6  G). 

The  union  of  all  intervals  in  IG,  considered  as  subsets  of  S,  defines  a  closed 
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subset  D  of  S,  since  IG  is  a  finite  collection  of  closed  sets  (10].  The  operation 


of  directed  rounding  will  now  be  defined  in  ID,  the  set  of  all  intervals  in  S  having 
endpoints  in  D. 

Definition  5.1.  For  x  E  d,  the  upward  rounding  operator  7  to  G  is  defined  by 

(5.2)  Vx  »  min  {b  |  b  2  x,  b  €  G)  , 
and  the  downward  rounding  operator  4  to  G  by 

(5.3)  4x  ■  max  {a(a£x,a£G}. 

For  X  •  (x  ,  x]  €  id,  the  directed  rounding  operator  □  to  IG  is  defined  by 

(5.4)  Qt  *  0(x  ,  x]  -  l Ax  ,  Vx]  ,  X  •  (JC  ,  x)  6  ID. 

It  follows  immediately  from  this  definition  that  □  is  an  inclusion  monotone  in¬ 
terval  operator.  Furthermore,  if  <J>  maps  D  into  itself  and  $  is  an  interval  extension 
of  4,  then  □$  will  be  an  interval  extension  of  $  which  maps  IG  into  IG.  Thus,  for 
actual  computation,  it  may  be  assumed  that  the  interval  extension  ♦  of  $  in  (2.3)  has 
been  constructed  to  map  IG  into  itself.  This  means  that  the  transformed  intervals 
«(Yn)  will  be  exactly  representable  in  terms  of  elements  of  G  for  Yn  £  IG.  The  theory 
of  interval  iteration  given  above  applies  to  operators  of  this  type  without  modifica¬ 
tion.  Furthermore,  under  the  following  reasonable  assumption,  the  entire  interval 
iteration  (2.3)  may  be  carried  out  exactly,  using  only  elements  of  the  grid  G. 

Assumption  5.1  (Intersection  Property).  If  X,Z  £  IG,  then  XnZ  »  0,  the  empty 
set,  or  XnZ  £  IG. 

Thus,  if  IG  has  the  intersection  property,  and  the  interval  extension  $  of  0 
has  been  constructed  to  map  IG  into  itself,  then  the  selection  of  an  initial  inter¬ 
val  Yg  £  IG  will  assure  that  the  intervals  Y^ ,  Y^ ,  ...  generated  by  the  interval  it¬ 
eration  (2.3)  also  belong  to  IG.  Furthermore,  this  interval  iteration  on  IG  will 
always  converge  or  diverge  in  a  finite  number  of  steps.  To  see  this,  let  G#X  - 
G# (x  ,  x]  denote  the  number  of  grid  points  (elements  of  G)  contained  in  the  interval 
X  €  IG.  Then,  the  following  result  holds. 
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Theorem  5.1.  On  1G,  the  interval  iteration  (2.3)  will  converge  or  diverge  in 
at  most  G# (Yq)  steps. 

Proof.  If  the  interval  iteration  (2.3)  has  not  terminated  in  n  steps  by  divergence 
or  convergence,  then  it  will  have  generated  distinct  intervals  YQ,  Y^ ,  . ..,  Y^  €  IG, 
including  the  initial  interval  YQ.  Since  (2.4)  holds,  one  has 

(5.5)  G#(Y  )  S  G#(Y  ,)  -  1  S  G#(Y  ,)  -2  i  . . .  S  G#  (Y  )  -  n. 

n  n-i  rw  0 

Thus,  the  maximum  length  of  a  sequence  of  distinct  nested  intervals  in  IG  is  v  = 

G#  (Yq)  -  1,  since  each  interval  in  IG  must  contain  at  least  one  point  of  G.  If  the 
interval  iteration  has  not  terminated  by  the  vth  step,  then 

(5.6)  G*lYu>  "  1 

by  (5.5)  ,  which  implies  that  Y ^  is  a  degenerate  interval.  Now,  one  has  either  Y^  C 
4(Y^),  in  which  case  Y^^  ”  Y^  (convergence  in  v  steps),  or  Y^n*(Y^)  =  0  (divergence 
in  v  +  1  “  G#  (Yq)  steps)  .  QED 

Thus,  in  actual  computation,  convergence  or  divergence  of  an  interval  iteration 
is  an  observable  event  in  principle,  since  one  works  on  a  grid  of  machine  numbers. 

Of  course,  G#  (YQ) ,  although  finite,  could  be  prohit Itively  large;  however,  termination 
of  interval  iteration  is  usually  observed  in  far  fewer  steps.  The  construction  cf 
the  interval  operator  4  is  crucial  to  the  success  of  interval  iteration  [3] ,  but  de¬ 
pends  heavily  on  the  nature  of  S,  $,  and  the  grid  G  available. 

6.  Solution  of  equations.  In  many  applications,  S  is  a  linear  space,  and  the  problem 
of  interest  is  to  find  a  solution  x*  of  the  equation 

(6.1)  f(x)  =  0. 

Thi3  may  be  transformed  into  a  fixed  point  problem  (1.1)  in  many  ways  by  the  intro¬ 
duction  of  suitable  iteration  operators  4-  For  example,  one  may  take 

(6.2)  *(x)  -  x  -  Yf (x) 

to  define  4,  where  Y  is  an  invertible  linear  operator  in  S.  Given  an  interval  ex¬ 
tension  F  of  f,  the  corresponding  interval  extension  4  of  $  is 
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(6.3) 


4(X)  =  X  -  YF (X) 


using  interval  arithmetic  [4],  [6].  If  S  is  a  Banach  space,  and  f  has  a  Fr£chet 
derivative  f',  then 

(6.4)  *' (x)  =  I  -  Yf* (x)  , 

where  I  denotes  the  identity  operator,  and  a  more  accurate  interval  extension  of  4 
can  be  constructed  on  the  basis  of  its  mean-value  form  [1] .  Let  F‘  be  an  interval 
extension  of  f‘,  then  the  corresponding  interval  extension  4*  of  4'  is  given  by 

(6.5)  4'  (X)  =  I  -  YF'  (X)  , 
and  the  mean-value  form  of  (6.2)  is,  for  y  £  X, 

(6.6)  4(X)  -  y  -  Yf  (y)  +  {I  -  YF'(X))(X  -  y)  , 

and  this  interval  extension  4  of  4  is  called  the  Krawczyk  iteration  operator.  It 
arose  from  consideration  of  am  interval  version  of  Newton's  method  (2) ,  and  has  many 
useful  properties  [3],  (51,  [9].  In  actual  practice,  computation  would  be  done  with 
a  rounded  version  of  (6.6).  Suppose  that  F,  F‘  are  interval  extensions  in  IG  of  f, 
f ' ,  then 

(6.7)  4{X)  -  Qly  -  YF (y)  +  {I  -  YF’(X)}(X  -  y)  } 

will  have  values  in  13  for  y  €  G,  X  £  IG.  An  even  more  rounded  interval  extension 
of  (6.6)  to  IG  is 

(6.B)  4(X)  =  □(□(y  -  []YF(y)}  +  0{D{I  -  OYF’  (X)  } -D(X  -  y)}}, 

which  models  a  realistic  computational  interval  operator. 
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